A Study of P-wave Heavy Meson Interactions in A Chiral Quark Model 
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The analytical forms of the interaction potentials between one S-wave and one P-wave heavy 
mesons as well as the potentials between two P-wave heavy mesons are deduced based on a chiral 
quark model. Our results explicitly show the attractive property between two heavy mesons. Con- 
sequently, a series of possible molecular states are obtained. It is expected that our study might 
shed some light on the popular discussions of the newly observed XYZ states. 
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The newly discovered XYZ states, such as new charmonium-like states of Y(3872) [l|-|4j, Y(3940) 
Y74050) ± X (4140) X(4160jjll], Y(4250) ± Y(4260) [HiL^, Y(4350) 0, Y(4360) 



H3, X(4430) ± 0-0], X(4660) 0,[ljand botomnium-like states of Z b (10610) ± 0, Z b (10650) 
21| . are of great interests, since they cannot be simply considered as normal quarkonium of cc or bb, 
especially for the charged ones. For these XYZ states, many different explanations have been taken into 
account, including tetraquark, molecular states, hybrid charmonia, cusps, and threshold effects. 

In our previous works 22J, l23| , the possible two S-wave heavy meson molecular explanations for the 
new resonances of Y(3872), Y(3940), ^(lOGlO) 1 * 1 and Zf,(10650) ± are obtained. In this paper we extend 
our study of the S-wave heavy meson interactions to the P-wave heavy meson interactions, and we try to 
explain the structures of the other XYZ states, particularly for ^(4250)*, Y(4350), Y(4360), Y(4260), 
Y(4430}±. As far as ^(4250)* is concerned, it was explained as a I G (J PC ) — + (l ) D* D\ molecule by 
Close [24J, however Nielsen et al. 25l427j| considered it as a 1 + (1~) DD\ molecule with QCD sum rules, 
while Ding 28( found that the DDi and DqD* molecules couldn't be the explanation for Y(4250) ± , but 



for X(4260). For the latter resonance. Close 24| and Nielsen et al. 2514271] both agreed with Ding's 
conclusion. Moreover, Nielsen et al. j25l427l ] concluded that Y(4430) ± could be regarded as a tetraquark 
or a D*Di molecule. Close [3 and Ding et al. [3] confirmed this molecular explanation, however, 
the calculation of Liu et al. [30( disfavored the D*D\ molecular explanation using the resonating group 
method in our chiral quark model. 

In this paper, our chiral quark model will be employed to study the heavy meson interactions. The 
model has been explained, in detail, in Refs. 31-38j]. In this chiral SU(3) quark model, onc-gluon-cxchange 
(OGE), confinement potential, scalar meson exchange as well as pseudoscalar meson exchange are taken 
into account. Moreover, this chiral SU(3) quark model was developed into the extended chiral SU(3) 
quark model by introducing the vector meson exchanges. It has been proved that both models are quite 
successful in reproducing the spectra of the baryon ground states, the binding energy of deuteron, the 
nuclcon-nucleon (NN), kaon-nucleon (KN) scattering phase shifts, and the hyperon-nucleon (YN) cross 
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sections 13 1H35I . During the past years, by using the Resonating Group Method (RGM), Liu 
and Wang [38| have studied the interactions and structures of the heavy-quark systems in both the chiral 
SU(3) quark model and the extended one. 
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Different from the RGM method, we have deduced a naly tical effective interaction potentials between 
the two S-wave mesons in our chiral quark models 22, 23, 3^|. In this paper, we try to go further to 
deduce the analytical effective interaction potentials between one S-wave and one P-wave heavy mesons 
as well as between two P-wave heavy mesons. Then, we'll systematically calculate the possible bound 
states associated with Dq , D\ , D\ and their bottom partners B\ , B\ , by solving the Schrodingcr equation 
with our analytical potentials, and try to explain the structures for some newly observed XYZ states, 
particularly for y>(4160), X(4250) ± , X(4350), X(4356), X(4260), and X(4430) ± . We expect our approach 
could give a more accurate description of the short-range heavy meson-meson interactions than the RGM 
does, especially comparing with Ref. (jol |. 

The paper is organized as follows. In section [TTl the framework of our chiral quark models is briefly 
introduced, and the analytical forms of the effective interaction potentials between one S-wave meson 
and one P-wave meson (S-P) and between two P-wave heavy mesons (P-P) in our chiral quark models 
are given. The bound state solutions are shown and discussed in Sec. IIIII Finally, a short summary is 
given in Sec. [IVJ 



II. THE CHIRAL QUARK MODEL 



The framework of our models has been discussed extensively in the literature 
mentioned in our previous papers 
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31-38, 40-4 
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As 



|, the internal kinetic energies and the internal interactions 
of each meson are not necessary to be calculated. The Hamiltonian of the relative motion, in our approach, 
takes the form of 



H = T, 



rel 



V, 



eff, 



(1) 



where T re i is the kinetic energy operator of the relative motion between the two mesons, and V e /f is the 
effective interaction potential derived from the quark-quark (quark-antiquark) interaction between two 
mesons by integrating out the internal coordinates £i and £2 of the two mesons: 



ij 3 



(2) 



while </?i(£i) and (^2(^2) are the intrinsic wavefunctions of the two mesons. If a meson is in S-wave, its 
wave function is taken as one-Gaussian form: 

¥>(0 = ( — ) e a 5 , 

and if in P-wave, it is 



(3) 



Vlm(0 



2%/2 5/4 _ 

d7T4 



'Y lm (0,<p). 



(4) 



Here, fi is the reduced mass of the two quarks inside each meson cluster and u is the harmonic-oscillator 
frequency of the meson intrinsic wavefunction. V(fij) in eq. ([2]) represents the interactions between the 
i-th light quark or antiquark in the first meson and the j-th light quark or antiquark in another. 

Because there is no color-interrelated interaction between the two color-singlet clusters, OGE interac- 
tion and confinement potential between two mesons don't exist, we only consider the meson-exchange 
interactions between the two meson clusters. In our chiral SU(3) quark model, we have 
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V(n 3 ) =£^(%) + £^(r%-) 



(5) 



a=0 



a=0 
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and in the extended chiral SU(3) quark model 



E^'(^)+E^-(f <3 -)+E^-(^). 



(6) 



V a "(fij) and V Ta (rij) are the scalar meson exchange and pseudoscalar meson exchange interactions, 
respectively. V pa (f f ij) indicates the vector meson exchange interactions. For quark-quark (antiquark- 
antiquark) interaction, the explicit forms of V a (7% ) , V^ a (fij) and V Pa (fij) have been described, in 
detail, in Refs. 
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31 



38 



40j-|43| as 



V^inj) = -C(g ch ,m aa ,A)X 1 (m aa ,A,r ij )(^X^), 

777. 

V^inj) = C(g ch ,m^ a ,A) — *° X 2 {m Va , A,r i3 ) fa ■ aj) (A? A?) , 



(7) 

(8) 



^ a (%) = C(5chv,rn Pa , A) 



X 1 (m Pa ,A,r ij ) + - 1 + — — J - + p / 

Qmitrij \ g chv M N g l h M N 



with 



x X2(m Pa ,A,rij) fa ■ aj) 



C(g ch ,m,A) 
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Xx(m,A,nj) = Yfanj) - —Y(Arij), 



X 2 (m,A,r ij ) 
Y(x) 



Y{mr i:j ) 
1 _„ 



r(Ar, 



— e 



(9) 

(10) 

(11) 
(12) 
(13) 



where A a is the Gell-Mann matrix in flavor space, and A is the cutoff mass indicating the chiral symmetry 
breaking scale. In eqs. (7-13), rrij and rrij are the masses of the i-ih light quark or antiquark in the first 



meson and the j'-th light quark or antiquark in another, respectively, while m 



and m Pa in eqs. (0 



[HI HI) are the masses of the scalar, pseudoscalar and vector nonets, respectively. Mm in eq. © is a mass 
scale usually taken as the mass of nucleon 33j. g c h, in the above eqs., is the coupling constants for the 
scalar and pseudoscalar nonets, while g c h v and f c hv are the coupling constants for the vector coupling 
and tensor coupling of vector nonets, respectively. For the quark-antiquark interactions, the G-parity of 
the exchanged mesons should also be taken into account. 

In the two heavy meson systems, we don't consider the one-meson exchange interactions between two 
heavy quarks or between one heavy quark and one light quark, because these interactions are beyond our 



SU(3) models. By using the method in Refs. 
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33 ] and integrating out the internal coordinates 



of two mesons as described by eq. ([2]), we can get the total analytical effective interaction potentials 
between two S-wave heavy mesons. Analogy to the two S-wave systems, we can go further to deduce 
the effective interaction potentials between one P-wave heavy meson and one S-wave heavy meson and 
between two P-wave Heavy mesons. These effective interaction potentials are listed as follows. 



Vejf (r) = E + E v & (*) + E v & («)> 
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with 



V qq a {R) = -Gr a C(g ch , mrTa ,A)X lq q(mr a ,A,R)(\%\1 
V£{R) = G„ a C(g ch ,m na) A) 



\2m q m q 



X 2qq {m« a ,A 1 R) (a q ■ a q ) (AgA|) , 



(14) 
(15) 



V q P i(R) 



G Pa C(g chv ,m Pa ,A) 
xX 2qq (m Pa ,A 1 R) (a q ■ a q ) 



n r 



1 + 



fchv m q + m q Jl hv m q m q 
gchv M N g 2 chv M% 



(16) 



Here, m q and m q are masses of the light quark and antiquark, respectively, G aa . 7Ta . Pa is the G-parity of 
the exchanged meson, and R is the relative coordinate between the two different mesons, namely, the 
relative coordinate between the centers-of-mass coordinates of the two mesons. Moreover in eqs. (14-16) 



X lqq (m,A,R) = Y qq (mR) - ^Y qq {AR), 
X 2qq (m,A,R) = Y qq (mR)-(-) Y qq (AR) 



where the Y qq (mR) is the modified Yukawa term. 

For the two S-wave heavy meson interactions, as mentioned in 22, 
in eqs. (|17I18|) reads 
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(17) 
(18) 

, the modified Yukawa term 



Y qq (mR) 



1 



2mR 



e -mR\ 1 _ er j 



P(R- 



2/3' 



1 - erf 



+ §)]}}• d9) 



For the one S-wave heavy meson and one P-wave heavy meson interactions, they receive the contributions 
from direct and exchange terms. The modified Yukawa term of the direct term reads 



Y qq {mR) 



1 



9 ? 



3mR I \ Acu(mQ + ^q) 2 ^ q Q 2 

'{l-erj[^{R+^)]}\-^Re-"> : - 



rnR I 



(20) 



and the one of the exchange term, associated with the charge-parity C in flavor wave functions (see eqs. 
(|27H28p in section III), reads 



Y qq (mR) 



3mR 



vhq d m q + mQ m q + tuq d 



m q + nig m q + ttiq da 4m,g tuq db 



)-i e tr<j e - mfl {l-er/ 
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Moreover, for the two P-wave heavy meson interaction, the modified Yukawa term is 



e 4jr 



2 2 



V4cj(toq + m q ) 2 iJ, q Q 2/ \4uj(mQ + rnqYuqQ 2/ 
-erf[- y/fi{R - } - e™*{l - er/ + g)] } 

TO?; 



x < e 



(m + m q ) 2 fi q Q w 2 (to q + mg) 2 ii q Q 
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-i?e 



In the above equations, 



Q 



(m Q + m q ) 2 n q Q (m Q + to 5 ) 2 /j 9 q 



(22) 



4a' 



where 



TOg + TOQ 



TOg + TOQ 



'to 9 + TOq\ 2 /to^ + to q n 2 

6 = MgQ^ ( ) + M?Q^ ' 



Q 



m Q 



m q + niQ rn Q rriq + toq m c 

c = /-V3 W — — ~ — 7~ZT~ ~ ^qQ u ' 



mQ rriq + m Q 



iriQ m q + mQ 



(23) 



(24) 

(25) 
(26) 



and m,Q and ijiq the are masses of the heavy quark and antiquark, respectively, and n q Q 



m q mQ 



In this work, we take our model parameters, which have been determined in our previous works, 



31-38, 40-42)]. The up/down quark mass m q is fitted as the nucleon mass and taken as 

2.621 is fixed 



as follows 

Mjv/3 ^313 MeV. The coupling constant for the scalar and pseudoscalar chiral fields g c h 
by the relation of 



An 



9 9NNn m u 

25 4tt M%'' 



with .9 2 v rAr 7r /47T = 13.67 determined from experiments. 

In our extended chiral SU(3) quark model, the vector coupling constant g c ] lv and tensor coupling 
constant f c h v in eqs. (|9l [T6|) arc fitted by the mass difference between N and A, under the condition that 
the strength of the OGE is taken to be almost zero. When the tensor coupling is neglected, g c h v = 2.351 
and f c hv = 0; and when the tensor coupling is considered, g c hv = 1-973 and f c hv = 1-315. The harmonic- 
oscillator frequency u> (being as 1/(to u o 2 ) = l/(m c b 2 ) where b u is fitted by the N — N scattering phase 
shifts) is taken as 2.522/m _1 in the chiral SU(3) quark model and 3.113/771" 1 in the extended chiral 
SU(3) quark model, respectively. In our calculation, the masses of the mesons are taken from the PDG 
44j . except for the a meson, which does not have a well-defined value. Here m a is obtained by fitting 
the binding energy of deuteron 33j. It is m a = 595 MeV in our chiral SU(3) quark model, 535 MeV 
for neglecting tensor coupling and 547 MeV for considering tensor coupling in our extended chiral SU(3) 



G 



quark model. The cutoff mass A means the chiral symmetry breaking scale and is taken as 1100 MeV as 
a convention. In addition, we find that the final results are not sensitive to the variation of the heavy 
quark masses, and we take m c = 1430 MeV 45[ and mj, = 4720 MeV 46 1 as the two typical values. 



III. RESULTS AND DISCUSSIONS 

To make sure the two- heavy meson systems have definite quantum numbers of isospin / and C-parity 
C, we follow the definitions of the flavor wavefunctions in Rcfs. 
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3G. 47 



For the hidden-charm 



systems constituted by two different mesons ID and ID, the flavor wavefunctions are 



r ^(S)'+S) + cS+S' ) 

Y 

1 = 1: < — (S)'-S° + cD-2)' ) 

r , - 

-[(ID °®° - ID +33") + c(S°ID - ID+D ")], 
1 = 0: i[(lD'°2) + ID'+ST) + c(£>°I)' + ©+a>' _ )]. 



(27) 



(28) 



Here, c=\ for C=+ and c=-l for C=— . As for ID meson-anti ID meson systems, the flavor functions are 



ID+S 
7 = 1: I 2TID 

— =(2)°2) - Q + Q~ 
v2 

i = 0: -Uid :d + S+S)-). 

V2 

For the two 25 systems, the similar expressions can be obtained. 



(29) 
(30) 



Two D masons 



1. S-P: DD 1} (DDI, DDI, D*D 1 ) and (D'*Dq, D* Dl) 



For the two pscudoscalar meson system of DD±, the spin factor < a q ■ er^ > in cqs. (Tl5l ITS)) is 0. 
Therefore, the spin-correlated interactions don't exist. Its interaction is associated with isospin I and 
C-parity C . We find two bound states, and their binding energies arc shown in Table [T] In the case of 
I = C = +, the I G (J PC ) = 0+(l-+) DDi bound state has a mass of 4253-4285 MeV. In the case of 
I = C = — , the total interaction potential is shown in Fig. [TJ and the 0~(1~~) DD\ bound state has 
a mass of 4264-4285 MeV, Its mass and quantum umbers match the X(4260). So in our approach, the 
0~(1~~) DDi molecule might be an explanation to the state of X(A260). 

For the systems of DDq, DD\ and D*Di, because they share the same interaction potential with the 
same isospin and C-parity, we put them together in discussion. We find the total interaction potentials 
in all cases are attractive, due to the strong attraction provided by a exchange. In the 1=0 C=+ case, 
the spin factor < a q ■ aq > in cqs. (|15[ I16[) is 1, and the attractive interaction potential, as shown in Fig. 
[21 is especially strong, because 7r and a' exchanges also provide strong attraction. 



r(fm) 



FIG. 1: The total interaction potential of 1=0 and C=— DDi system. The solid, dashed and dotted lines represent 
the results obtained from the chiral SU(3) quark model, the extended chiral SU(3) quark model neglecting and 
including the tensor coupling of vector field, respectively. 

o 
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l - 10 ° 
> 

-150 
-200 
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FIG. 2: The total interaction potential of 1=0 and C=+ DDq, DD2 and D* D\ systems. The solid, dashed and 
dotted lines represent the same meaning as in Fig. [Tj 

Not only in this 1=0 C=+ case but also in the case of 1=0 C=- can we find the bound states after 
solving the Schrodinger equation. We list the obtained binding energies in Table U and one may find the 
binding energies, for the systems of DDq, DD% and D*D\, are almost the same, and the small differences 
are due to the small differences among their reduced masses. 

We find the + (0 h ) DDq molecule has a mass of 4154-4181 MeV, and it could be a plausible expla- 
nation for the resonance of X(4160). The mass of 0~(0 ) DDq molecule, in our calculation, is about 
4160-4183 MeV. Moreover, we get the 0+(2 _+ ) DD* 2 molecule with the mass of 4297-4325 MeV, and 
it could be an explanation for X(4350). The 0"(2 — ) DD\ molecule has a mass of 4303-4326 MeV, 
and the + (J~ + ) D*D X or 0"(J"") D*D X molecules has a mass of 4400-4429 MeV or 4406-4430 MeV, 
respectively. However, no 1=1 D*D\ molecule could be found in our model, because 7r and a exchanges 
are both repulsive. Therefore, our result doesn't favor the D*D\ molecular explanation of the isovector 
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TABLE I: Binding energies (MeV) of the bound states of S-P D meson systems. I, II and III refer to the chiral 
SU(3) quark model, the extended chiral SU(3) quark model neglecting and considering tensor coupling of vector 
field, respectively. 





I G (J PC ) 




I 


II 


III 


DDi 






0.4 


32.4 


22.6 








0.3 


21.0 


12.0 


DDn 


+ (0" + ) 




1.7 


29.3 


19.6 




0~(0~~) 




0.4 


23.2 


15 


DDI 


0+(2-+) 




1.8 


29.7 


20.0 




0"(2~) 




0.5 


23.6 


15.4 


D*D X 


0+(J"+) 




2.2 


31.2 


21.2 




o-(j-) 




0.7 


24.7 


16.3 


D*D* 


o+(i-+)(s= 


-2) 


21.8 


30.7 


30.3 




o-(i~ )(S= 


-2) 


5.9 


12.6 


12.1 


D*D* 2 


+ (J~+)(S= 


=2) 


22.3 


31.3 


30.9 




o-(j— )(S= 


=2) 


6.1 


12.9 


12.4 



X(4430) ± . This result is consistent with the conclusion of Liu 30(, and contrary to the results of Close 



2j and Ding [2 



For the two different vector meson systems, such as D*Dq and D*D\, their interaction potentials are 
simultaneously associated with the isospin /, the C-parity C and the total spin S. In the 1=0, C=+, 
5=2 case, the spin factor < a q -Uq >, in eqs. ([T5l[TfI| . is 2, which is different from that of 1=0 C=+ of the 
DDq, DD\ and D*D 1 systems. We find that the D*Dq and systems share the same interaction 

potential shown in Fig. [31 and correspondingly we get the bound states with the binding energies listed 
in the Table HI In this case the effects of the vector meson exchanges are small in the extended chiral 
SU(3) quark model. In the case of 1=0, C=-, S=2, the bound states also exist. 

In our calculation, the mass of 0+(l~+)(S=2) D*D^ molecule is of 4297-4306MeV and we find the 
(T(l — )(S=2) D*Dl molecule has a mass of 4315-4322MeV, and it might be explained as X(4360) 
No bound states could be found in the 1=1 cases in our approach. Our result agrees with Ding 



28 



that X(4250) ± might not be explained as a isovector D*D$ molecule. Our results also disfavor the 
D*Dq molecular explanation of X(4260) 25l427|. Since we ignore the spin-orbital coupling interaction in 
our calculation, the + (J~+)(S=2) D*D* 2 molecule with mass of 4441-4452MeV and the 0"(J — )(S=2) 
D*D 2 molecule with mass of 4459-4466MeV are degenerate states for the total angular momentum J. 



2. P-P: DiDi, (DiDq, DiD* 2 ) and (D* D* Q , D* 2 D* 2 , D* D* 2 ) 

For the system of D\D\, its interaction is only associated with isospin /. It could form a + (J ++ ) 
weakly bound state with binding energy listed in Table [TTJ This bound state has a mass of 4824-4842 
MeV in our approach. The interaction potential is shown in Fig. 01 
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-200 - 



1 2 

r(fm) 

FIG. 3: The total interaction potential of 1=0 C=+ S=2 D* Dq and D* D% systems. The solid, dashed and dotted 
lines represent the same meaning as in Fig. [T] 
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FIG. 4: The total interaction potential of 7=0 -Di-Di system. The solid, dashed and dotted lines represent the 
same meaning as in Fig. Q] 

For the D\Dq and T)\D\ systems with the same quantum numbers, they share the same interaction 
potential. Also in the 1=0 C=+ case the interaction potential becomes strong enough to form bound 
states. This interaction potential is shown in Fig. [5j We see that the interaction potential in the extended 
chiral SU(3) quark model is weaker than that in the chiral SU(3) quark model, because the vector meson 
exchanges, in this case, provide repulsive interaction, however the range of the total potential is longer 
than that of the chiral SU(3) quark model. As a result, the vector meson exchanges don't obviously 
change the binding energies as listed in Table [TT1 In our approach, the + (l ++ ) DiD^ and + (J ++ ) 
D x Dl molecules have the masses of 4719-4727 MeV and 4863-4870 MeV, respectively. 

For the two vector systems of DqDq, D^D^ and DqD^, they share the same interaction corresponding 
to the same isospin I and total spin S, and the possible bound states are also listed in Table [TT1 In the 
case of 1=0 S=2 DqD'q, D\D\ and DqD^ systems, the spin factor < <r q ■ <7q > in cqs. ([Pol [TB|) is 1, 
the same as that of the 1=0 C=+ D\Dq and D\D\ systems. Here the interaction potential is shown 
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FIG. 5: The total interaction potential of 7=0 and C=+ D 1 Dq and D\D 2 systems and 1=0 5=2 DqDq, D2D2 
and D0D2 systems. The solid, dashed and dotted lines represent the same meaning as in Fig. [T] 

TABLE II: Binding energies (MeV) of the bound states of P-P D meson systems. I, II and III refer to the same 
meaning as in Table [I] 





I G 


(J PC ) 




I 


II 


III 




0+ 


(J ++ ) 




0.4 


18.2 


12.0 


D1D0 


0+ 


(1++) 




12.2 


19.5 


12.4 


D 1 D* 2 


0+ 


(J ++ ) 




12.5 


19.8 


12.4 


D* D* 


0+((T 


-+)(S= 


2) 


12.1 


19.3 


12.4 




0"(0" 


~)(S= 


1) 




16.7 


11.4 




0+((T 


-+)(S= 


0) 




15.6 


11.3 


D$Dl 


0+(J 


f+ )(S= 


=2) 


12.7 


20.0 


12.9 




o-(j- 


^")(S= 


=1) 




17.5 


12.1 




o+(j- 


f+ )(S= 


=0) 




16.3 


11.9 


D* D* 2 


0+(2^ 


-+)(S= 


2) 


12.4 


19.7 


12.6 




0+(2" 


-+)(S= 


1) 




17.1 


11.7 




0+(2" 


" + )(S= 


0) 




16.0 


11.6 



in Fig. [5l which is strong enough to form the bound states in our models. Our results tell the possible 
0+(0 ++ ) L>o-Dq molecule with the mass of 4617-4624 MeV, the 0+(J++) D\T)* 2 molecule with the mass 
of 4904-4911 MeV and the possible 0+(0++) D*D* molecule has a mass of 4760-4768 MeV. 

For the systems of 7=0 5=0,1 DqDq, D^D^ and D^D^, we get the bound states in our extended 
chiral SU(3) quark model, but no bound states in our chiral SU(3) quark model. This is because, in 
our extended chiral SU(3) quark model, vector meson exchanges provide additional strong attractive 
interaction. Even though we can't currently draw a definite conclusion whether or not the molecular 
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bound states could exist, those states provide a test for our two models. 

As a short summary, we have studied twelve sets of the heavy meson interactions associated with 
P-wave D mesons. They are, S-P: DD X , {DDI, DD* 2 , D*D X ) and (D*D^, D* D* 2 ); and P-P: D X D U 
(DiDq, DiDl) and (DqDq, D 2 D 2 , DqD 2 ). We've found eighteen possible bound states, and they are 
all isoscalar. Only the four possible molecules of 0+((T + ) DD%, 0"(1 — ) D*D* Q , 0"(1~) DD U and 
+ (2~ + ) DD* could be explained as V(4160), X(4360), A(4260), V(4350), respectively Other bound 
states don't match any observed XYZ states. 

B. Two B masons 

Considering the resemblance between B mesons and D mesons, the corresponding binding systems 
should have the similar properties, while the P-wave B meson bound states should have larger binding 
energies because of their larger reduced masses. 

1. S-P: BB\, (BB2, B\B*) and B* B% 

As for the BB\ system, similar to the DD\ system, it becomes bound in the I = C = ± cases and 
the obtained binding energies are listed in the Table Hill In the I = C = + case, the + (l h ) BB\ 
molecule has a mass of 10932-10980 McV, and in the I = C = — case, the _ (1~~) BBi molecule has 
a mass of 10955-10993 MeV, 

For BB 2 and B\B* , the systems of one S-wave B meson and one P-wave B meson, they share the same 
interaction potential corresponding to the same isospin and C-parity. In the 1=0 case the interaction 
potential is strong enough to bind BB 2 and B\B* . For example, in the 1=0 and C=+ case, the interaction 
potential is shown in Fig. HI and the obtained binding energies are listed in Table HID We find that the 
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FIG. 6: The total interaction potential of 1=0 C=+ BB'2 and B\B* systems. The solid, dashed and dotted lines 
represent the same meaning as in Fig. [1] 

+ (2~ + ) BB 2 and 0"(2 — ) BB* 2 molecules have the mass of 10962-11009 MeV and 10977-11018 MeV, 
respectively. In our calculation the + (J h ) B\B* molecule, without a definite total angular momentum 
J, has a mass of 10970-11017 MeV, and the mass of _ (J~) B X B* molecule is about 10985-11026 MeV. 
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TABLE III: Binding energies (MeV) of the bound states of S-P B meson systems. I, II and III refer to the same 
meaning as that in Table [I] 







I G {J PC ) 


I 


II 


III 


BBi 




o+(i-+) 


20.0 


69.0 


56.1 






o-(i~) 


7.7 


46.1 


34.5 


BB$ 




+ (2" + ) 


17.8 


64.6 


51.2 






0~(2~~) 


9.2 


50.1 


38.6 


B\B* 




o+(j-+) 


17.8 


64.7 


51.3 






Q-(J-) 


9.2 


50.2 


38.6 


B*B* 2 


H 


-(J-+)(S=2) 


66.8 


97.8 


91.7 




o- 


-(J-+)(S=1) 


2.2 


9.4 


7.7 




0" 


-(J-)(S=2) 


24.9 


40.8 


40.1 



For the _B*£?2 system, the total interaction potential is associated with isospin /, C-parity C=+ and 
total spin 5. The possible bound states are listed in the Table ITO! The interaction potential between B*B\ 
and in the case of 1=0 C=+ 5=2 is shown in Fig. [7] This interaction potential is strong enough to bind 




FIG. 7: The total interaction potential of 1=0 C=+ S=2 B*B% system. The solid, dashed and dotted lines 
represent the same meaning as in Fig. [1] 

B*Bl system, and the 0+(J~+)(S=2) B* molecule has a mass of 10974-11005 MeV in our calculation. 
In the case of 1=0 C=+ 5=1, the mass of 0+(J _+ )(S=l) B*B* molecule is about 11063-11070 MeV, 
while in the case of 1=0 C=- 5=2, the 0~(J— )(S=2) B*B\ molecule has a mass of 11061-11047 MeV. 
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2. P-P: Si Si, S1B2 and B* 2 Bl 

For the two pseudoscalar meson system of B±Bi, the total spin S is 0, like DiD ll the total interaction 
potential is only associated with the isospin /. Like the D\D\ system, B\B\ forms a bound state in the 
1=0 case in our models as shown in Table HVl Our + (J ++ ) B\B\ molecule has a mass of 11408-11436 
MeV. 

For the B\B 2 . like the D\D 2 system, the interaction potential shown in Fig. [51 in the case of 1=0 and 
C=+, is strong to bind the B\B 2 . The obtained binding energies in this case are listed in the Table HVl 




r(fm) 

FIG. 8: The total interaction potential of 7=0 and C=+ SiS| system and 1=0 5=2 D2D2 system. The solid, 
dashed and dotted lines represent the same meaning as in Fig. [T] 

We find that the 0+(J++) B X B* 2 molecule has a mass of 11435-11444 MeV. 

TABLE IV: Binding energies (MeV) of the bound states of P-P B meson systems. I, II and III refer to the same 
meaning as in Table [I] 







i g (j pc : 




I 


II 


III 


Si Si 




o+(j++; 




7.1 


39.5 


30.5 


B1B2 




o+(j++; 




31.1 


33.2 


24.3 






o-(j+-; 






20.2 


14.2 


SJ B2 


0^ 


(J++)(S= 




31.2 


32.5 


23.7 




0" 


(J+-)(S 


=1) 




36.8 


30.4 




0^ 


(j++)(s 


=0) 




40.1 


35.7 



As for the system of B^B^, the interaction potential is associated with isospin / and the total spin S. 
Similar to the D^D^ system, we find a bound state in the 1=0 5=2 case where the interaction potential 
is just the same as the 1=0 C=+ B\B 2 system shown in Fig. [5J The binding energies are also listed in 
the Table[lVl In our approach, the 0+(J++)(S=2) molecule has a mass of 11454-11465 MeV. 

To summarize this subsection, we have studied seven sets of heavy meson interactions associated with 
P-wave B mesons. They are S-P: BB^, (B*B U BBi), and B*B* 2 and P-P: B X B* 2 , B* 2 B* 2 , and B X B X . 
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We've found twelve possible bound states and they are all isoscalar. Unfortunately, none of these bound 
states could match any newly observed XYZ states. 



IV. SUMMARY 



In this work, we have employed our SU(3) chiral quark model to systematically studied the heavy 
meson interactions associated with the P-wave. The parameters of our models have already fixed in our 
previous calculations from the spectra of baryon ground states, the nucleon-nucleon phase shifts and 
the deuteron binding energy. We have gotten some possible molecular states. For the systems of P-wave 
B-mesons, unfortunately, those molecules don't match to any of the newly observed XYZ hadronic states. 
For P-wave D meson interactions, we find that the X(4350) could be explained as the + (2 '") DD2 
molecule, the A(4360) could be explained as the 0~(1 ) D*Dq molecule, and the A(4160) could be 
explained as the + (l h ) DD^ molecule. However, the X(4430) ± could not be explained as a isovector 
DiD* charged molecule. Moreover, neither X(4250) ± nor A(4260) could be explained as the D*Dq 
molecule, but A(4260) might be explained as a _ (1~~) DD\ molecule. 

Finally, we find that we can get bound states in our extended chiral SU(3) quark model in some cases 
but not in the chiral SU(3) quark model. This is because the vector meson exchanges in the extended 
chiral SU(3) quark model provide additional attraction. Future experimental results about these special 
cases would give a discriminates of our two models. 
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